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Real-time dynamics of charge density and lattice displacements is studied during photoin- 
duced ionic-to-neutral phase transitions by using a one-dimensional extended Peierls-Hubbard 
model with alternating potentials for the one-dimensional mixed-stack charge-transfer complex, 
TTF-CA. The time-dependent Schrodinger equation and the classical equation of motion are 
solved for the electronic and lattice parts, respectively. We show how neutral domains grow 
in the ionic background. As the photoexcitation becomes intense, more neutral domains are 
created. Above threshold intensity, the neutral phase is finally achieved. After the photoexci- 
tation, ionic domains with wrong polarization also appear. They quickly reduce the averaged 
staggered lattice displacement, compared with the averaged ionicity. As the degree of initial 
lattice disorder increases, more solitons appear between these ionic domains with different 
polarizations, which obstruct the growth of neutral domains and slow down the transition. 
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1. Introduction 

Dynamics in non-equilibrium phase transitions has 
. been challenging and studied from different viewpoints. 
In cooperative processes, roles played by microscopic de- 
grees of freedom are interesting. When photoirradia- 
' tion triggers a phase transition, induced real-time dy- 
, namics gives important information on their roles. We 
' here focus on photoinduced phase transitions in quasi- 
, one-dimensional mixed-stack charge-transfer (CT) com- 
' plexes, which have been studied extensively both exper- 
, imentally 1-8 ) and theoretically. 9 ~ n ) 

Quasi-one-dimensional mixed-stack CT complexes 
consist of donor (D) and acceptor (A) molecules, which 
■ are alternately stacked to form columns. They exhibit 
a neutral-ionic phase transition in equilibrium by ap- 
plying hydrostatic pressure or by lowering the temper- 
ature. 12 '' In these complexes, electrons are itinerant and 
an electron-lattice coupling plays a key role in causing 
, the phase transition. At the phase transition charge is 
transferred between D and A molecules, accompanied by 
a distortion forming (D +s A~ s ) dimers. The two phases 
are distinguished from each other by the ionicity [p). It 
is a measure of the charge transfer and defined later so 
that the neutral (N) phase has p < 0.5, while the ionic 
(I) phase has p > 0.5. 

Tetrathiafulvalene-p-chloranil (TTF-CA) is one of 
such CT complexes and shows a photoinduced phase 
transition. 1-8 ) In the TTF-CA complex, the TTF 
molecule is a donor, and the CA molecule is an accep- 
tor. 13 ) At ambient pressure, it is ionic below Tni = 81 K 
and neutral above Tni- In the ionic phase, two structural 
phases exist. When the inversion symmetry is lost, it 
becomes ferroelectric. 14 ) It is realized at ambient pres- 
sure in equilibrium. When the polarizations D+^A - " 5 
and A~ S D +S are disordered to recover the symmetry, 



it becomes paraelectric. It is realized under high pres- 
gure i5) or j ug ^. a ft er photoexcitation. 8 ) The pressure- 
temperature phase diagram is similar to typical solid- 
liquid-gas phase diagrams. 

Recently, the time-resolved reflectivity has been mea- 
sured with improved resolution, which reflects the evo- 
lution of the spatially averaged ionicity. 5 ~ 7 ) Time- 
resolved second harmonic generation (SHG) has also 
been measured to study the evolution of the degree 
of the inversion-symmetry breaking after photoexcita- 
tions. 8 ) The SHG intensity decreases rapidly after the 
photoirradiation of the ferroelectric ionic phase, but the 
ionicity decreases much slower. This information is im- 
portant to reveal non-equilibrium processes in the pho- 
toinduced phase transition. Because the reflectivity and 
the SHG intensity reflect the averaged properties, we 
do not know creation or annihilation processes of neu- 
tral and ionic domains. Space-resolved properties of the 
coupled charge-lattice system are important to clarify 
the non-equilibrium processes during the photoinduced 
phase transition. 

Theoretically, Koshino and Ogawa have studied pho- 
toinduced structural dynamics in a one-dimensional sys- 
tem composed of classical lattice displacements and lo- 
calized electrons occupying one of the two levels per 
site. 16 ) The dynamics is described by nucleation of a 
kink-antikink pair and motion of the kinks. The strength 
and the range of the elastic-lattice coupling have been 
shown to decide whether a global change is realized or 
not. Meanwhile, Huai et al. have studied the sta- 
bility of a neutral domain in the ionic background in 
a quasi-one-dimensional itinerant-electron-lattice model 
for TTF-CA. 10 ) The origin of threshold excitation in- 
tensity for the transition is discussed, based on the adi- 
abatic potential as a function of the domain size. Its 
dynamics was not studied there. Only few theoretical 
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studies have so far been devoted to dynamics in pho- 
toinduced phase transitions of itinerant-electron-lattice 
systems. 11, 17 - 18 ) Our previous and preliminary studies 
have shown only that the motion of the ionicity has three 
characteristic time scales, which originate from the elec- 
tronic motion, the lattice motion, and the neutral-ionic 
domain-wall motion, during the ionic-to-neutral phase 
transition. 17, 18 ) 

In this paper, we study the dynamics of coupled 
charge density and lattice displacements during the pho- 
toinduced ionic-to-neutral phase transition in a one- 
dimensional extended Peierls-Hubbard model with alter- 
nating potentials for TTF-CA. We show how neutral do- 
mains grow or shrink in the ionic background. The dy- 
namics is described by domain walls, which exist between 
ionic and neutral domains or between ionic domains with 
different polarizations. In §2, the model and the method 
are described. In §3, results are presented. Finally in §4, 
the conclusion is given. 

2. Extended Peierls-Hubbard Model with Al- 
ternating Potentials 

The highest occupied molecular orbital (HOMO) of 
the D molecule and the lowest unoccupied molecular or- 
bital (LUMO) of the A molecule are necessary to de- 
scribe the electronic state. In the limit of vanishing 
transfer integrals to - ► 0, the HOMO is doubly occu- 
pied in the neutral phase, while both of the HOMO and 
the LUMO are singly occupied in the ionic phase. Dif- 
ferent types of electron-lattice couplings have so far been 
employed to explain the dimerization. 19-22 ) Here we as- 
sume that the Coulomb interaction strength is modified 
by the lattice displacement. 22 ) We use a one-dimensional 
extended Peierls-Hubbard model with alternating poten- 
tials at half filling, which is a single-chain analogue of the 
model used in rcf. 10, 
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where, c\ a (ci i<T ) is the creation (annihilation) opera- 
tor of a 7r-electron with spin a at site I, ni t „ — cj a ci t(T , 
n i = n L] + n i,l; u i is the dimensionless lattice displace- 
ment of the Zth molecule along the chain from its equidis- 
tant position, and y\ = ui+i — it;. The distance be- 
tween the Ith and (I + l)th molecules is then given by 



ri = ro(l + ui + i — ui), where r is the averaged distance 
between the neighboring molecules along the chain. The 
D and A molecules are located at the odd and even sites, 
respectively. The nearest-neighbor repulsion strength 
between the ith and (I + l)th sites Vi depends on the 
bond length yi, Vi — V + foyf, where V is for the reg- 
ular lattice, and (3 2 is the quadratic coefficient. The lin- 
ear coefficient is neglected for simplicity. The parameter 
to denotes the nearest-neighbor transfer integral, U the 
on-site repulsion strength, and d the level difference be- 
tween the HOMO and the LUMO in the neutral phase 
with to — 0. The elastic energy is expanded up to the 
fourth order: the parameters fci and ki are the linear and 
nonlinear elastic constants. The mass of the Ith molecule 
is denoted by mi. The number of sites are denoted by 
N. 

2.1 Hartree-Fock approximation 

We use the unrestricted Hartree-Fock (HF) approxi- 
mation, which replaces the electronic part (2.2) by 

N 
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where h^ a is the Hartree term, hf is the Fock term, 
o~i,even = 1 for even I, and Si yeven — otherwise. We itera- 
tively solve the eigenvalue equation to obtain the ground 
state that is self-consistent with the one-body densities 
(c; a ci t(T ) and (c\ a ci + i^) and the lattice displacements 
yi, which satisfy the Hcllman-Feynman theorem. 

We then add random numbers to the initial yi and ui 
values according to the Boltzmann distribution at a ficti- 
tious temperature T. They roughly correspond to ther- 
mal lattice fluctuations just before the photoexcitation. 
The degree of initial lattice disorder T used here is much 
smaller than the energy barrier between the neutral and 
ionic phases. With thus modified lattice displacements, 
the HF Hamiltonian is diagonalized once again by impos- 
ing the self-consistency on the one-body densities only. 
Photoexcitations are introduced by changing the occu- 
pancy of the highest occupied and lowest unoccupied HF 
orbitals. 
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2.2 Equations of motion 

For the electronic part, we solve the time-dependent 
Schrddinger equation, 



\^, a (t)) =Texp 



|^(0)>, (2.8) 



where T represents the time-ordering operator. For tj 
j x At < t < tj + At, it is rewritten as 



\ip„, a (tj + At)) =Texp 



tj+At 



dt'H* F (t>) 
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The exponential operator is decomposed, 23 ) in such a 
way that the decomposition is accurate to the order of 
(At) 2 . 24 ) 

For the lattice part, we solve the classical equation of 
motion, 



mi- 



dt 2 
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3.1 Adiabatic potentials 

Figure 1 shows the adiabatic potentials in the ionic 
phase (dotted line) and in the neutral phase (solid line), 
as a function of the staggered lattice displacement y st , 
which is defined as yi = (— 1) y s t- The molecules are 
equidistant for y st — 0, and dimcrized for \y st \ =/= 0. The 
neutral phase has a local minimum at y s t = 0, while the 
ionic phase has two minima at about y s t = ±0.06. In the 
ionic phase, the donor and acceptor molecules have pos- 
itive and negative charge, respectively. In the dimerized 
•lattice, the polarization vectors are uniform so that this 
phase is ferroelectric. With the parameters above, the 
ferroelectric ionic phase is stable, and the neutral phase 
is metastable. The initial state in the later calculations 
is always in the ferroelectric ionic phase. The optical 
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where the expectation value is taken with respect to the 
wave function (2.8). We employ the leap frog method, 
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which are accurate to the order of (At) 2 . The computa- 
tion time is shorter than the similarly accurate, velocity 
Vcrlet method. 25 ^ The lattice velocity at time t can be 
defined as 



vi(t) 



Vl(t+^)+Vl(t-^) 



(2.14) 



During the numerical calculations we have confirmed 
that the total energy is always conserved. 

3. Results and Discussions 

We use A=100, f =0.17 cV, [7=1.528 eV, ^=0.604 
eV, d=2.716 cV, /3 2 =8.54 eV, fc!=4.86 eV, fc 2 =3400 eV, 
and the bare phonon energy oj = (l/r )(2fci/m r ) 1 / 2 =0.0192 
eV, and impose the periodic boundary condition. Here 
the reduced mass m r is defined as m r — tod"1a/('^d + 
toa) with the mass of the donor molecule mo and the 
mass of the acceptor molecule m\- With these parame- 
ters, the dimerized ionic phase is stable and the neutral 
phase is metastable. These parameters are almost the 
same as in ref. 10, but the bare phonon energy w used 
here is about four times higher than the optical phonon 
energy of the TTF-CA complex in the neutral phase. 
Note that ui here is in the unit of the DA distance in the 
regular neutral phase. The averaged ionicity is defined 
asp=l + l/ATE^ 1 (-l)'(^). 



-0.68 



Fig. 1. Adiabatic potentials in the neutral phase (solid line) and 
in the ionic phase (dashed line), as a function of the stag- 
gered lattice displacement y s t- The parameters are to =0.17 
eV, U =1.528 eV, V =0.604 eV, d =2.716 cV, /3 2 =8.54 eV, 
fci =4.86 cV, and k 2 =3400 cV. 

phonon energy in the ionic phase is 2.4 times larger than 
that in the neutral phase. 

3. 2 Emergence of a few I-I solitons after weak photoex- 
citations. 

We show here what happens in the ionic phase af- 
ter weak photoexcitations. When two electrons are ex- 
cited and the degree of initial lattice disorder is very 
small, T—0. 00017 eV, the ionicity and the staggered 
lattice displacements evolve as shown in Fig. 2. The 
horizontal axis gives the elapsing time t multiplied by 
u. The solid line in Fig. 2 (a) shows the averaged 
ionicity p. After the photoexcitation, it first decreases 
from its initial value p = 0.92 to p = 0.76 and it soon 
comes back to the initial value by the time about 10/w. 
The space and time dependence of the boundaries be- 
tween the ionic and neutral domains, where pi = 0.5, is 
shown in Fig. 2 (b). Here the ionicity pi is defined as 

Pi = ^^{-(ni-i) + 2{m) - (ni +1 )) + 1. The vertical 
axis gives the site index I. It is clearly seen that neutral 
domains are created immediately after the photoexcita- 
tion and they soon disappear by the time about 10/w, 
which coincides with the behavior of p. The dotted line in 
Fig. 2 (a) shows the magnitude of the averaged staggered 
lattice displacement defined as y st = l/iV| X^ =1 (— l) 1 Vil- 
li would roughly be proportional to the SHG intensity. 
Its initial value is y st = 0.06. At first, y Bt quickly drops 
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Fig. 2. (a) Averaged ionicity p (solid line) and averaged staggered 
lattice displacement y s t (dashed line), as a function of time t af- 
ter the photoexcitation, which is multiplied by the bare phonon 
energy ui. (b) Space and time dependence of the boundaries be- 
tween the neutral and ionic domains (p; = 0.5). (c) Correlation 
between the staggered lattice displacement y st ; and the ionic- 
ity Pi, as a function of the site index I and the elapsing time 
t multiplied by u>. Inset (I) of (c) shows the area indicated by 
the box on the enlarged scale. The parameters arc to =0.17 
eV, U =1.528 eV, V =0.604 cV, d =2.716 cV, fa =8.54 eV, 
fci =4.86 cV, k 2 =3400 cV, u> =0.0192 cV, and T =0.00017 cV. 
In the 100-site chain, 2 electrons are excited. 



to y st — 0.035 much before the recovery of p at the time 
about 10/oj. Then, y st slightly increases on average to 
y st — 0.045 by the time about 20/uj. Finally, y st grad- 
ually decays after that. The bars in Fig. 2 (c) show the 
correlation between the staggered lattice displacement 



{-Vl-i + ZjJi - 2/l+i) an d the ion- 



dchned as y st i - 4 
icity pi, as a function of the site index I (vertical axis) 
and the elapsing time t multiplied by 10 (horizontal axis) . 
As shown in Fig. 3, each bar has two components. The 
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Fig. 3. Long and slanted bar (a) on distorted ionic sites, and short 
and horizontal bar (b) on undistorted neutral sites, in Fig. 2 (c). 
The vertical component shows y a t 1, while the horizontal one pi. 

vertical component shows y s t 1 , while the horizontal one 
pi. Then, distorted ionic sites are described by long and 
slanted bars [Fig. 3 (a)], while undistorted neutral sites 
by short and horizontal bars [Fig. 3 (b)]. In ionic do- 
mains with (D + A~) [(A~D+)] dimers, y st 1 is positive 
(negative), so that the bar points up (down) [Fig. 3 (a), 
solid (dotted) bar] . These domains are denoted by I and 
I domains, respectively, hereafter. It is clearly seen that, 
after the photoexcitation, neutral domains with almost 
equidistant molecules arc quickly created at first. This is 
why y s t decays in the corresponding period in Fig. 2 (a). 
Then, small I domains appear after the neutral domains 
disappear [Fig. 2 (c)] , so that y st is not recovered to the 
initial value any more. The boundary between the I and 
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Fig. 4. I-I solitons between I and I domains with different polar- 
izations, which are denoted by the directions of the bars. D and 
A arc donor and acceptor molecules, respectively. 



I domains is regarded as a I-I soliton, which looks like 
Fig. 4 in Fig. 2 (c). Finally, the number of I-I solitons 
gradually increases. The excitation energy is partially 
absorbed into these solitons. In this subsection, we show 
the creation of small neutral domains by weak photoex- 
citations and their decay into I-I solitons. 
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3.3 Photoinduced ionic-to -neutral phase transition 

In this subsection, we show how neutral domains grow 
in the ionic background after strong photoexcitations. 
When four electrons are excited and the degree of initial 
lattice disorder is still very small, T=0. 00017 eV, the 
ionicity and the staggered lattice displacements evolve 
as shown in Fig. 5. The seed used to produce random 
numbers is the same as that in the previous subsection 
and will be used later throughout the paper. The mean- 
ings of the lines are also the same as before. After the 
photoexcitation, the averaged ionicity p first shows a dip 
and gradually decreases on average to p = 0.27 by the 
time about 45/w [Fig. 5 (a), solid line]. The magnitude 
of the averaged staggered lattice displacement y st de- 
creases on average from its initial value y st = 0.06 to 
y st ~ by the time about 30 /u [Fig. 5 (a), dotted line]. 
Two neutral domains are created by the photoexcitation, 
but the smaller one immediately disappears [Fig. 5 (b)]. 
Meanwhile, the larger neutral domain expands and anni- 
hilates the initial ionic domain by the time about 45/w. 
It is clearly seen in Fig. 5 (c) that, when the smaller neu- 
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Fig. 5. (a) Averaged ionicity p (solid line) and averaged staggered 
lattice displacement y s t (dashed line), as a function of time t af- 
ter the photoexcitation, which is multiplied by the bare phonon 
energy u>. (b) Space and time dependence of the boundaries be- 
tween the neutral and ionic domains (p; = 0.5). (c) Correlation 
between the staggered lattice displacement y st i and the ionicity 
pi, as a function of I and uit. Insets (I) and (II) of (c) show the 
areas indicated by the boxes on the enlarged scale. The param- 
eters are the same as in Fig. 2. In the 100-site chain, 4 electrons 
are excited. 

tral domain disappears, I and I domains appear instead. 



Because of the appearance of the I domains, the aver- 
aged staggered lattice displacement y s t decreases much 
faster than the averaged ionicity p. The expansion of the 
larger neutral domain is shown in the insets (I) and (II) 
of Fig. 5 (c) on the enlarged scale. In the inset (I), the 
expansion is smooth because both y st i and pi change 
gradually. In the inset (II), the expansion is stepwise 
due to I-I solitons. The I-T solitons near the neutral- 
Neutral domain 
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Fig. 6. I-I solitons near the neutral-ionic domain boundary ob- 
struct the growth of the neutral domain. 



ionic domain boundary obstruct the growth of the neu- 
tral domain, as shown in Fig. 6. When the small but 
finite staggered displacements around the neutral-ionic 
domain boundary become parallel to those of the ionic 
domain, the neutral domain rapidly expands. In this 
subsection, we show the obstruction of the growth of the 
neutral domain by I-I solitons. 

3.4 Case of many I-I solitons 

In this subsection, we show a case of many I-I soli- 
tons created by the photoexcitation before the ionic-to- 
neutral phase transition. The number of I-T solitons sen- 
sitively depends on of the degree of initial lattice disorder 
T: it increases with T. When four electrons are excited 
and the degree of initial lattice disorder is not very small, 
T"=0.0017 cV, the ionicity and the staggered displace- 
ments evolve as shown in Fig. 7. The averaged ionicity 
decreases from its initial value p ~ 0.878 to p ~ 0.266 
by the time about 130/cj [Fig. 7 (a)]. Here three char- 
acteristic time scales are observed. The averaged ionic- 
ity p rapidly oscillates. We calculated the charge-charge 
and current-current correlation functions in the random 
phase approximation (RPA). 26 ' Their spectra have peaks 
at about 6 x 2ttlo and 5 x 2iru>, respectively, in both the 
neutral and the ionic phases at T = 0. The correspond- 
ing periods are then about 0.16/cj and 0.2/o>, respec- 
tively. A period of about 0.16/w is indicated by "EL" in 
the inset of Fig. 7 (a). The clearly seen beat is due to 
the oscillation of period about 0.16/w from the charge- 
charge correlation and that of period about 2/u) from 
the current-current correlation. The slow oscillation of p 
indicated by "PH" in Fig. 7 (a) is correlated with that 
of y s t- Thus, the observed periods ranging between 3/w 
and 6/ui are due to the lattice oscillations. It is noted 
that the period of the noninteracting lattice oscillation is 
2tt/lo. Recall that the optical phonon energy in the ionic 
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Fig. 7. (a) Averaged ionicity p (solid line) and averaged staggered 
lattice displacement y s t (dashed line), as a function of time t af- 
ter the photoexcitation, which is multiplied by the bare phonon 
energy ui. (b) Space and time dependence of the boundaries be- 
tween the neutral and ionic domains (pi = 0.5). (c) Correlation 
between the staggered lattice displacement y at ; and the ionicity 
pi, as a function of I and cut. Inset (I) of (c) shows the area 
indicated by the box on the enlarged scale. The degree of initial 
lattice disorder is T =0.0017 eV. The other parameters are the 
same as in Fig. 2. In the 100-site chain, 4 electrons are excited. 



phase is 2.4 times larger than that in the neutral phase. 
Thus, the observed range of the period distribution is 
reasonable. The slowest time scale is due to the collec- 
tive motion of domain walls and indicated by "DW" in 
Fig. 7 (a). A growth of the neutral domain is observed 
during the period between and 130/w, accompanied 
with the shrinking ionic background and with the oscil- 
lating domain boundaries [Fig. 7 (b)]. The decrease of p 
in Fig. 7 (a) is due to the growth of the neutral domain in 
Fig. 7 (b). The averaged staggered lattice displacement 
y st decreases much quicker than p [Fig. 7 (a)]. After 
the photoexcitation, many I domains appear [Fig. 7 (c)]. 
It suggests that the SHG intensity becomes weak very 
quickly, compared with the decrease of the ionicity, even 
in the one-dimensional model. Many I-I solitons are ob- 
served in the inset (I) of Fig. 7 (c), and they indeed 
slow down the growth of the neutral domain. This sit- 
uation is already illustrated in Fig. 6. The I-I solitons 
are obstacles to the expansion of the neutral domain, 
and they sensitively affect the collective motion of the 
neutral-ionic domain walls. 

Fourier transforms of some time windows of the ionic- 
ity p in Fig. 7 (a) are shown in Fig. 8 (a). Here, the dot- 
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Fig. 8. (a) Fourier power spectra for the averaged ionicity p in 
Fig. 7 (a). The dotted line is obtained during the I-to-N tran- 
sition (0/uj < t < 102/u;), and the solid line after the transi- 
tion (143/oj < t < 245/w). (b) Charge-charge (solid line) and 
current-current (dotted line) excitation spectra in the RPA for 
the static neutral state. 

ted line is obtained from < t < 102/ 'u during the I-to-N 
transition, and the solid line from 143/w < t < 245/u 
after the transition. During the transition, slow com- 
ponents are dominant and broad reflecting the complex 
motion of the neutral-ionic domain walls coupled with I- 
I solitons. After the transition, the broad spectrum has 
been attenuated, and rather sharp peaks appear instead. 
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For the static neutral state, linear excitation spectra are 
calculated in the RPA and compared with these Fourier 
spectra. The solid and dotted lines in Fig. 8 (b) show 
those of the charge-charge correlation and of the current- 
current correlation, respectively. Here, the charge den- 
sity difference between the odd and even sites (wave vec- 
tor 7r) and the total current density (wave vector 0) are 
respectively used in the RPA calculations. The peak 
at energy about 6.2 x 2iru) in Fig. 8 (a) corresponds to 
the charge-charge correlation. The position of the lower- 
energy peak in Fig. 8 (a) depends on the number of ex- 
cited electrons, and it is shifted to the lower-energy side 
with the increasing number (not shown) . Finally, the po- 
sition of the lower-energy peak comes close to the peak 
of the current-current excitation spectrum in Fig. 8 (b). 
This is interpreted as follows. After the phase transition, 
the total electronic energy is still higher than the total 
energy of the metastable, static and undistorted neutral 
state. The excess energy creates many electron-hole pairs 
in the neutral phase. The number of electron-hole pairs 
increases with the excess energy, and thus with the num- 
ber of excited electrons. When the number of electron- 
hole pairs is large enough, the pairs are overlapped with 
each other. These overlaps overcome the somehow dis- 
ordered lattice displacements so that they enhance the 
collective nature of electron-hole pairs, i.e., the excitonic 
nature. The exitonic effect reduces the energy of charge 
transfer and slightly slows down the motion of the ion- 
icity. It is now clear that the rapid oscillation of the 
ionicity is composed of two modes, one of which is of the 
excitonic origin. 

3. 5 Dependence on the number of excited electrons 

We show here the dependence on the number of ex- 
cited electrons. When six electrons are excited and the 
degree of initial lattice disorder is not very small again, 
T=0.0017 eV, the ionicity and the staggered displace- 
ments evolve as shown in Fig. 9. The averaged ion- 
icity p decreases faster [Fig. 9 (a)] than that after the 
four-electron excitation [Fig. 7 (a)]. Now there are many 
neutral- ionic domain boundaries [Fig. 9 (b)]. The neu- 
tral domains grow much faster than those after the four- 
electron excitation, so that the ionic-to-neutral phase 
transition is completed by the much shorter time about 
27 /u. Some I-T solitons are seen in Fig. 9 (c). The neu- 
tral domain grows so rapidly that its motion is insensitive 
to the I-I solitons. 

The dependence of the evolution of the averaged ion- 
icity p on the number of excited electrons with the 
degree of initial lattice disorder, T = 0.0017 eV and 
T = 0.00017 eV, is shown in Figs. 10 (a) and 10 (b), 
respectively. The dashed, dotted, and solid lines show p 
after two (denoted by Ex2 in the figures), four (Ex4), and 
six (Ex6) electrons are excited, respectively, in the 100- 
site chain. Here, the time required for the completion of 
the photoinduced phase transition is denoted by T n Tm, 
where n describes the number of excited electrons and 
to distinguishes the degree of initial lattice disorder. As 
more electrons are excited, the transition is completed 
earlier. It means that the stronger photoexcitations in- 
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Fig. 9. (a) Averaged ionicity p (solid line) and averaged staggered 
lattice displacement y s t (dashed line), as a function of time t af- 
ter the photoexcitations, which is multiplied by the bare phonon 
energy ui. (b) Space and time dependence of the boundaries be- 
tween the neutral and ionic domains (p; = 0.5). (c) Correlation 
between the staggered lattice displacement y st ; and the ionicity 
pi , as a function of / and u)t. The parameters are the same as in 
Fig. 7. In the 100-site chain, 6 electrons are excited. 
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Fig. 10. Dependence of the evolution of the averaged ionicity p on 
the number of excited electrons, (a) with T =0.0017 eV, and 
(b) with T =0.00017 eV. The dashed, dotted, and solid lines 
show p after two (denoted by Ex2), four (Ex4), and six (Ex6) 
electrons are excited, respectively, (c) shows the dependence on 
the degree of initial lattice disorder. The dashed, dotted, and 
solid lines show p after 4 electrons are excited with T =0.00017 
eV, T =0.00085 eV, and T =0.0017 eV, respectively. The other 
parameters arc the same as in Figs. 2, 5, 7, and 9. 



duce the faster ionic-to-neutral transitions. The depen- 
dence on the degree of initial lattice disorder T is shown 
in Fig. 10 (c). The dashed, dotted, and solid lines show 
p after four electrons are excited with T = 0.00017 eV, 
T = 0.00085 cV, and_ T = 0.0017 eV, respectively. As 
T increases, more I-I solitons are produced. Because 
the I-I solitons obstruct the growth of neutral domains, 
they slow down the transition. Consequently the tran- 
sition period r„Tm becomes longer as the initial lattice 
variables are more strongly disordered. Disorder itself 
is necessary to nucleate neutral domains. At T = 0, 
neutral domains can be created by numerical errors, but 
they do not appear immediately, so that the transition 
is not always completed quickly. We remark that, in 
real systems, the free-energy difference and the barrier 
height between the ionic and neutral phases also depend 
on the temperature. As long as the initial randomness 
in the lattice displacements and velocities are concerned, 
higher temperatures are unfavorable to reduce the tran- 
sition period. However, the other effects remarked here 
would work on the contrary. 

4. Conclusions 

We have investigated how neutral domains grow or 
shrink in the ionic background after photoexicitations 
through their dependence on the excitation intensity and 
on the degree of lattice disorder in the initial state, using 
the one-dimensional extended Peierls-Hubbard model. 
As the photoexcitation becomes intense, more neutral 
domains are initially created and their average size in- 
creases. Above threshold intensity, it eventually becomes 
neutral throughout the system. After the photoexci- 
tation of the ionic state, neutral domains, I domains, 
and thus I-I solitons are created. The I domains (i.e., 
with the opposite polarization to the initial one) reduce 
the averaged staggered lattice displacement y st - Conse- 
quently, the averaged staggered lattice displacement y st 
decays much faster than the averaged ionicity p, in gen- 
eral. This result is consistent with the experimentally 
observed, quick decay of the SHG intensity, 8 ) compared 
with the decrease of the ionicity, although the calculation 
is performed in the one-dimensional system. As the ini- 
tial lattice variables are more strongly disordered, more 
I-I solitons appear. The IT solitons obstruct the growth 
of neutral domains. Therefore, disorder slows down the 
phase transition. In short, the progress of the photoin- 
duced ionic-to-neutral phase transition is not simply like 
dominoes that are falling down, in the sense that the 
competition between I-I solitons and neutral-ionic do- 
main boundaries affect the growth of neutral domains 
and thus the transition itself. 
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